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L  INTRODUCTION 


The  electromagnetic  scattering  from  an  arbitrarily  shaped  shell  of  revolution  has  been 
investigated.  An  example  of  the  geometry  is  shown  in  Fig.  1 .  The  surface  is  parameterized  in 
the  cylindrical  coordinate  by  (p„(5),  (j).  with  0  <  5  <  /  and  0  <  cf)  <  2ti  ,  where 

S  o 

the  z-axis  is  the  axis  of  rotation  and  s  is  the  arc  length  parameter  of  the  generating  curve  (Fig. 
2).  Henceforth  the  subscript  “g”  will  denote  the  described  on  the  generating  curve.  Also 
shown  in  Fig.  2  are  the  unit  tangent  f  and  the  unit  normal  ^  ^  x  f  to  the  curve.  The 

parameterization  of  the  generating  curve  is  chosen  so  that  « *  is  the  outward  normal  of  the 
outer  surface  of  the  shell.  The  angle  0^,  measured  from  the  positive  f  axis  to  f ,  unlike  the 
polar  angle  in  a  spherical  coordinate  system,  can  assume  negative  values  or  positive  values 
greater  than  7t. 

The  inner  and  outer  surfaces  of  the  perfectly  conducting  shell  are  coated  with 
anisotropic  materials  which  can  be  different.  It  is  assumed  that  the  thickness  of  the  coating 
and  of  the  shell  are  infinitesimal.  The  coated  surfaces  are  also  assumed  to  satisfy  the 
impedance  boundary  condition  (IBC).  It  is  found  that,  with  some  special  impedance  matrices, 
the  backscattering  cross  section  along  the  axis  of  this  coated  shell  can  be  eliminated.  This 
finding  will  be  presented  in  this  report. 

For  convenience  in  both  theoretical  formulation  and  numerical  computation,  it  is 
desirable  to  consider  £  as  the  electric  field  intensity  vector  divided  by  the  intrinsic  impedance  =y[^Jz 
of  the  isotropic,  homogeneous  medium  within  which  the  shell  is  located.  Therefore,  electric 
field  intensity  jv  takes  the  unit  of  amperes  per  meter,  the  same  as  those  of  the  magnetic  field 
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intensity and  the  electric  and  the  magnetic  surface  current  distributions.  The  surface 
impedances,  Z*  on  the  outside  and  Z '  on  the  inside  of  the  shell,  are  also  normalized  with  the 
same  factor  to  become  dimensionless.  The  time  dependence  ^  will  be  assumed  and 
suppressed.  These  are  the  conventions  adopted  in  this  report. 

For  the  remainder  of  this  chapter,  the  major  steps  taken  to  prove  the  vanishing  of  on- 
axis  backscattering  cross  section  are  briefly  sketched  in  Sections  A  through  D.  The  sufficient 
condition  to  eliminate  this  cross  section  is  specified  at  the  end  of  Sec.  D.  In  Sec.  E,  the 
impedances  with  which  the  on-axis  backscattering  can  be  eliminate  are  given.  They  are  the 
main  contributions  of  this  work.  All  general  purpose  programs  [1,  2]  for  the  numerical 
computation  of  scattering  cross  sections  of  anisotropic  impedance  coated  objects  should  be 
checked  for  their  accuracy  against  selected  surfaces  of  revolution  coated  with  impedance 
matrices  satisfying  the  conditions  listed  in  that  section.  The  results  will  provide  an  indication 
of  the  error  bound  for  the  particular  algorithm. 


A.  SCATTERED  FIELD 

The  Stratton-  Chu  equations  [3]  give  the  radiation  in  a  homogeneous  and  isotropic 
region  outside  a  closed  surface  S  in  terms  of  the  tangential  components  of  the  £  and  ff  field 
intensities  on  S.  The  scattered  fields  £  and  at  a  point  r  outside  S  which  encloses  the 
shell  are  given  by: 

4TiE'"(r)  =  /cV  X  f  [nxE(rJ]G(r-rJda^+  ik^  f  [nxff(FJ]G(r-FJda^ 

s  J  s 

-  iV  f  [nxB(rJ]-V^G(F-rJda^ 

s 


J 


(1) 


ATzH'\f)  =  kv  X  r  [/}xB(rJ]G(r-rJda^-  ik^ 

s 

-f  iV  f  [nxE(rJ].Vfi(r-rJda^ 

^  s 

ik\f-fj 

where  n  is  the  unit  outward  normal  of  S,  ^  =  cov/ils,  and  G(r-Fj=— - .  Note  that  the 

k\f-f  I 

tangential  components  of  £  and  b  ^^e  the  equivalent  magnetic  and  electric  surface  current 
distributions  on  S  (rotated  by  90*^ .) 

In  eqs.  (1)  and  (2),  the  surface  S  can  be  shrunk  toward  the  shell.  We  denote  the  outer 
surface  of  the  shell  as  and  the  inner  surface  of  the  shell  as  5'".  Henceforth  the 
superscripts  “+”  and  will  denote  the  described  values  on  S  ^  and  S "  respectively.  The 
total  (incident  and  scattered)  equivalent  electric  currents  on  ^  ^  and  on  5  “  are; 


f  {nxEirJlGif-Qda^ 

(2) 


K"  =  n^xB^ 

=  n~  xB~  =  -n^xB 


(3) 


where  fi "  and  n  '  are  the  unit  outward  normals  on  S  *  and  S  ~  respectively,  g  ^  and  -  are 
the  total  magnetic  field  intensities.  Because  the  thickness  of  the  shell  is  considered 
infinitesimal,  n*  =  -fi  '  at  any  point  of  the  shell. 

Similarly,  the  equivalent  magnetic  currents  and  are: 


C  =  E"  X  n* 

L  =  E  X  n'  =  E  X  (-fi  ’■) 


(4) 


Define  .^((j)^,  5^)  as  the  sum  of  the  outside  and  inside  total  electric  surface  currents,  and 
similarly  ^  ^  as  the  sum  of  the  outside  and  inside  total  magnetic  surface  currents; 
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K  =  +  K~  (5) 

L  =  r  +  L'  (6) 


Since  the  thickness  of  the  coating  and  of  the  shell  are  assumed  to  be  infinitesimal,  the 
contribution  from  and  Si  (Fig.  3)  to  the  integrals  must  vanish  by  the  edge  condition  [5] 
as  the  distance  of  S  to  the  shell  becomes  zero.  Eqs.  (1)  and  (2)  are  reduced  to: 


and 


47r.£  ^‘^(r) 

=  -  kV 

x/' 

Jg  . 

r 

JQ 

^,sJG{f-rJp^ds^d(\>^ 

+  ik^ 

f'r 

JO  JO 

)G(f-rJpdsd(^^ 

-,Vj 

flr2n 

0  Jo 

-e- 

X 

> 

II 

IT 

JG  JO 

■0)  G(r-f  )p/5/(l)^ 

+  ik^ 

r  1  rlix 

io  Jo 

1  G(r-fJp^ds^d(J>^ 

-,V| 

'  1  r2n 

0  JO 

0  Vfi{r-r;)pdsd<^^ 

(7) 


(8) 


Note  that  the  scattered  fields  are  determined  completely  by  the  sum  currents  on  the  shell. 

In  the  far-field  region,  the  scattered  field  can  be  expressed  in  the  spherical 

coordinate  system  as: 

E'\r)  ~  pj  0[^sin0^cos0cos((j)-(j)  )+Z,cosesin((t)-(l)  ) 

-.^,cos0„sin0  +i:^cos(({)  -4>^)  -Z-^sin0gSin(({)  -(j)^)] 

+4>[^^cos(({)-(l)J-i:,sin0gSin(4)-4)J  (9) 

-Z,jSin0gCos0cos(<|)  -cj)  J  +Z-,cos0„sin0 

-Z^cos0sin(4)-({)^)] 

where  the  subscripts  cf)  and  t  denote  the  tangential  components  of  the  surface  currents  in  the  ^ 
and  f  directions  on  the  surface  of  the  shell  respectively.  Because  of  the  rotational  symmetry 
of  the  shell,  the  Fourier  expansion  can  be  utilized  to  solve  the  (j)  -dependence  of  this  problem. 
Define  the  Fourier  series  component  (5)  of  a  function  /  ((t),5)  by; 
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(10) 


/(<t>4  =  E  e“*/„  (s) 

n~-oo 


then  as  0  -  0,  eq.  (9)  is  simplified  and  transformed  to: 


E\f)  = 


[(K^,  ^iL^,)smQ^](x^iy) 


(11) 


The  backscattering  cross  section  is  given  by; 


E  ^ 

Or.  =  lim  - 


(12) 


For  an  incident  plane  wave  directed  along  the  axis  of  the  shell,  with  normalized  to 

unity,  Og  becomes: 


^5  = 


f  Po^ 

Jo 


~ikz^ 


4(^<^,_r^^<t,-i)-^(^,-i-^4,-i)sinej(f-/T)}|^ 


(13) 


B.  IMPEDANCE  BOUNDARY  CONDITION 

The  impedance  boundary  condition  links  the  tangential  components  of  the  total 
electric  field  to  those  of  the  total  magnetic  field  on  a  surface  through  its  surface  impedance 
which  is  a  function  of  the  electromagnetic  properties  of  the  material,  i.e.. 
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X  n^)  =  Z"  (n^  X  H^) 


(14) 


n-  X  (£* 


where  Z  ^  and  Z '  are  the  normalized  surface  on  5  ^  and  S '  respectively.  Equivalently, 

«=  X  L"=  Z=  R-  (15) 


For  anisotropic  materials,  it  is  more  convenient  to  utilize  matrix  notations.  In  eq.  (15), 
the  vectors  and  £-  can  be  considered  as  two-element  column  vectors,  with  each  of  the 
(})-  component  designated  as  element  1  and  the  t  -  component  as  element  2.  and  Z  “  are 
then  two-by-two  matrices.  The  cross  product  of  the  unit  vector  n  *  with  such  a  two-element 
column  vector  can  be  represented  with  the  following  matrix  multiplication; 


n^x  L"  =  V V 


(16) 


where  V  = 


0  -1 

1  0 


.  Note  that  V 


T  _ 


-V  and  V-  =  -I. 


In  terms  of  the  matrix  V,  the  impedance  boundary  condition,  eq.  (15),  can  be  written 


as: 


Z-K  ^-V  L 


(17) 


Since  the  scattered  field  is  determined  completely  by  the  sum  currents,  it  is  desirable 
to  write  the  impedance  boundary  condition  as  one  which  provides  the  difference  currents  in 
terms  of  the  sum  currents.  Define: 
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(18) 


Z  = 


-  Z-) 


and 


A  =  I  (Z-  -  Z  )  ,19) 

then,  assuming  that  Z  is  invertible,  eq.  (17)  can  be  transformed  into  the  following  two 
equations. 


-  K'  =  Z-^VL  -  Z'^A  K 


(20) 


L^-  L'  =  -V[Z  -  AZ’^A]  K  -  V[^Z-^Vi  (21) 


C.  EQUATIONS  FOR  SCATTERING  CURRENTS 

From  eq.  (3),  using  matrix  notation  and  noting  that  fi  x  H  =  n  x  H  on  the 

tan 

surface  of  the  shell,  one  can  deduce  that: 


K 


K'  =  -  {-V)H 


tan 


nC'  -C)  -  nsr  +  C')  =  y  (HZ’  *  hZ  *  2hZ) 


(22) 


therefore. 


+ 

■‘■‘tan 


2HZ'  =  -V{r  -K-) 


(23) 


Similarly, 
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(24) 


eZ:  -  ez  -  2EZ  -  V (r  -  L-) 

The  Fourier  coefficients  of  the  tangential  components  of  the  scattered  fields 
approaching  the  outer  and  inner  surfaces  of  the  shell  in  eqs.  (23)  and  (24)  can  be  obtained 
from  eqs.  (7)  and  (8)  and  written  as: 

^  =  -K  EZ)  +  Z  (25) 


Kn(^)  -  -K  -  K 


(26) 


The  elements  of  the  two-by-two  matrices  and  are  integrodifferential  operators  on  the 
elements  of  the  column  vectors  and  i^.  They  are  derived  and  given  later  in  eqs.  (85)  and 
(86). 

Combining  these  with  eqs  (20)  and  (21); 


K 


K 


VZ'^A 


AZ-^F 


-  ' 

K 

me 

h 

n 

lan,/7 

> 

=  2 

£ 

ginc 

n 

(27) 


where  unlike  the  products  of  M,,  iV„  with  Z’  Z’  >^atrices  involving  the  impedances  are 
constants  and  the  product  with  are  usual  matrix  multiplications.  Given  the  incident 

wave,  this  is  a  linear  integral  equation  of  the  second  kind  for  the  unknown  currents  K  ,  L  ■ 
It  is  assumed  that  unique  solutions  to  eq.  (27)  exist  for  all  n. 
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D.  ZERO  ON-AXIS  BACKS CATTERING  CONDITION 


Because  of  the  rotational  symmeti^,  a  wave  incident  along  the  axis  of  the  shell  (the 
z-axis)  can  be  assumed  to  be  linearly  polarized  in  the  y-direction  without  loss  of  generality. 
Let  the  fields  on  the  surface  of  the  shell  be: 

=  y  e  [  cos(j)(^  +  sin(})  (sinG^f  +  cosG^fi"")  ]  e 

to  to 

(28) 

’”‘^(4),5')  =  X  e  =  [-sincf)^  +  cos<f>  (sinG^  +  cosG^fi")  ]  e 


where  the  superscript  “wc”  denotes  the  incident  wave. 

The  Fourier  coefficients  of  the  tangential  components  of  the  incident  fields  are; 


-ikz„ 


t^n 


Kn  =  [S„-,  - 


le  ^ 


-ikz 

k:  = 


(29) 


Note  that  only  the  ^  =  ±1  terms  in  eq.  (29)  are  nonzero.  Therefore,  only  the  cases 
n  =  ±\  have  to  be  considered. 

Since 


pmc 
^tan,±  I 


1 

e  ^ 

1 

+7sinG 

sj 

tan,±l 

2 

sinG^ 

-ikz^ 


(30) 


The  incident  fields  satisfy  the  relation: 


11 


0 

-il 

1 

il 

0 

^inc 

;^tan,±l 

/ 

0 


(31) 


On  the  other  hand. 


0  -il 

M 

n 

n 

1 

0  -il 

il  0 

N 

n 

M 

n 

il  0 

(32) 

Therefore,  if 


Z-AZ-*A 


VZ'^A 


AZ'^V 


also  commutes  with 


0 


il 


-il 

0 


then 


Z-AZ'^A 

AZ-'F 

\ 

0  -il 

^*1 

Hi 

1 

> 

-VZ-^V 

Ir^ 

il  0 

0  -il 

II 

Hi 

-Hi' 

Z-AZ~^A  AZ’H 

Hi 

il  0 

11 

Hi 

Hi. 

+ 

1 

<I 

>• 

.4 

1- 

0 

-il 

} 

[ 

il 

0 

1 

^  taii,±l 

(33) 


The  existence  and  uniqueness  of  the  solution  to  eq.  (27)  implies  that 


1 


0  -il 
il  0 


=  0 


(34) 
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i.e.,  ^  ±  =  Q .  From  eq.  (13),  a^=  0.  There  will  be  no  backscattering. 


E.  IMPEDANCE  MATRICES 

It  has  been  deduced  that  a  sufficient  condition  to  eliminate  on-axis  backscattering  is: 


0  -ii\ 

Z-AZ'^A  AZ^^Z 

AZ'^Z 

0  -il 

il  0 

ZZ‘*A  -VZ-^V 

<1 

-ZZ'^Z 

il  0 

which  is  equivalent  to: 


=  -AZ''F 


(36) 


and 


Z+  FZ-iF=AZ-'A 


To  satisfy  these  two  conditions,  either 


Z"  =  Z'  = 


0 

+i 


±i 

0 


(37) 


(38) 


when  Z  ^  and  Z '  are  skew-symmetric,  or,  in  the  case  when  Z  ^  and  Z '  are  symmetric, 

detZ'"  =  detZ’  =  1  (39) 

det(Z"  -  Z')  =  0  (40) 
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Note  that  when  the  shell  is  a  closed  surface,  the  impedance  boundary  condition 
separates  the  interior  volume  from  the  solution  in  the  exterior.  Only  the  conditions  on  Z  ^  will 


apply  to  the  exterior  problem.  Therefore,  for  an  anisotropic  impedance  coated  body  of 

0  ±/ 

revolution,  the  on-axis  backscattering  cross  section  vanishes  either  -  or 

+7  0 

detZ"  =  1  with  a  symmetric  Z^ 
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n.  SCATTERmC  OF  A  COATED  SHELL  OF  REVOLUTION 


A.  RADIATION  FROM  EQUIVALENT  CURRENTS 

The  geometry  and  coordinate  system  for  an  arbitrarily  shaped  shell  of  revolution  are 
shown  in  Figures  1  and  2.  In  the  cylindrical  coordinates  (p,  (|),  z),  the  shell  is  described  by 
p  =  Pg(5),  z  -  zj^s)  with  0<5</  where  s  is  the  arc  length  parameter  of  the  generating  curve. 
Since  both  p  and  z  are  functions  of  5,  the  surface  of  the  shell  is  parameterized  in  the 
coordinates  cf)  and  s.  The  transformations  of  vectors  between  the  cylindrical  coordinates  and 
the  surface  coordinates  (4>,  s)  are  given  in  Appendix  A. 

The  Stratton-Chu  equations  give  the  radiation  in  a  homogeneous  and  isotropic  region 
outside  a  closed  surface  in  terms  of  the  £  and  ^  field  intensities  on  that  surface.  Using  the 
conventions  adopted  in  this  report,  the  radiated  fields  and  at  a  point  f  outside  S 
which  encloses  the  shell  of  revolution  can  be  written  as  [4] : 

=  ik'^  f  (n><ff(rJ)G(r-rJda^  +  k  f  («x£(T))  x  V^G(r-fJda^ 

(41) 

+  k  f(n  •  E(rJ)  Vfi(r-fjda^ 

s 

4TiH'Xr)  =  -  ik^  f  (nxE(rJ)G(r-Fjda^  +  k  f  x  Vj3{r-r^)da^ 

^  S  J  s 

,  -  («) 
*k  {n>  H(rJ) 

s 


where  k  -  coyjae  and  G{f-f^  =  — - ,  and  n  is  the  outward  normal  of  S. 

k\r-r^\ 
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Eqs.  (41)  and  (42)  can  be  converted  into  ones  involving  integrals  of  tangential 


components  of  the  fields  only.  Consider  the  last  term  in  eq.  (41): 


j  {ft  0  E) 


WG  da^  = 

o  o 


■  jin  «  E) 


VG  da 


=  -iVJ 

^  (fi 

s 

®  ^G  da^ 

II 

1 

< 

n  ® 

(V. 

X  ^G  da^ 

=  -,Vj 

r  n  < 

’  [V, 

X  (^  G)  - 

=  -,v/ 

■[V,x 

(H  G)]dv^ 

=  -/  V/ 

in 

s 

X  H) 

o  VJ3  da^ 

(43) 


Note  that  VG(f-f^  -  “V^G(r-f^).  The  second  term  in  eq.  (41)  can  be  written  as: 

j(n>‘E(0)  X  Vfi(F-rJ  da^  ^  ^  x  j(nxE(rj)  G{f^rJ  da^  (44, 

Consequently,  eq.  (41)  is  converted  to: 

ATiE'Xr)  =  A  V  X  r[n  X  E(rJ]G(r-rJda^+  iP  (  [n  ^  H{fJ\G{f-r^da 

S  •'  iS 


i  V  X  H(fJ\  0  Vfi{r-f^da^ 


(45) 


From  the  duality  principle  [6],  eq.  (42)  can  be  obtained  from  eq.  (41)  by  replacing  ^  with  g 
and  }j  with  _  £.  The  corresponding  expression  to  eq.  (45)  for  the  radiated  magnetic  field  can 
be  obtained  similarly: 

4-KH'\f)  =  ^  V  X  X  H{fJ\G{r-f^da^-  iP  j  \n  x  E{fJ\G{f-f^da^ 

+  /■  V  I  [n  X  £(f  )]  «  Vfi(r-rjda^ 
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In  eqs.  (45)  and  (46),  the  surface  S  consisting  of  the  four  non-overlapping  surfaces 
S^,  S^,  S^,  and  S2  can  be  shrunk  to  the  surface  of  the  shell  so  that  S'"  and  -  S~ 
(Fig.  3).  Since  the  thickness  of  the  coating  and  of  the  shell  are  assumed  to  be  infinitesimal,  the 
edge  conditions  [5]  guarantee  that  the  contributions  to  the  integrals  from  the  bottom  surface 
and  the  top  surface  must  vanish  in  this  limit  (Appendix  B).  Therefore,  the  integrals  in  eqs. 
(45)  and  (46)  need  to  be  carried  out  over  S and  S  ‘  only. 

The  total  (incident  and  scattered)  equivalent  electric  currents  ^  *  on  S"  and  ^  on  ' 
are  defined  to  be; 


K*  =  n^xH^ 
K'  =  fi-xH^  =  -n^xH 


(47) 


and  the  equivalent  magnetic  currents  on  ,5"  and  on  ‘  are  defined  as; 


C  =  E"  X 

L'=  E'  xn-  =  E'x(-n^) 


(48) 


With  each  integral  in  eqs.  (45)  and  (46)  written  as  a  sum  of  integrals  over  the  two  surface  S 
and  5",  we  can  substitute  the  surface  currents  of  eqs.  (47)  and  (48)  for  the  fields. 
Furthermore,  define  the  sum  electric  current  as; 

K  =  +  K~  (49) 


and  the  sum  magnetic  current  as; 
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L  =  L  +  L 


(50) 


eqs.  (45)  and  (46)  are  reduced  to; 


4tzE  '\r) 


kV  X  r  G(r-0 

V  0  J  Q 

iP  ^(4><,,0  G(f-fJ 

V  0  0 

zV  /■'/■"”  ^((i)^,5^)  •  V^G(f-f ) 
»'  0  0 


(51) 


4nH-\r)  =  kV  x  G(r-fJ 

*  ik-  f'E  m^.sj  C(r-rj 

*'0  *'0 

-/V  YY’'  r((l)^,5^)  ®  Vfi(f-Fj  pdsd4>^ 

d  0  J  0 


(52) 


where  f  =  (p„,(f),,^J  and  p„  =  pi5^),  =  zXsJ. 


B.  THE  FAR  FIELD 

In  the  far-field,  the  spherical  coordinate  system  (r,  0,  (}>)  is  more  convenient.  To  the 

lowest  order  in  Y ; 

r 


G  (r-rj 


ikr 

6  £  ^  /:r^cos0+  /:p^sm0cos(4) -(()^)  ] 

r 


(53) 


and 
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2,  0 

where  cos6  =  -  and  sin6  =  -t-.  Eq.  (51)  simplifies  to: 
r  r 

E"\r)^ik- — f  ds^f  V({)^pJ0[i:,sin0gCOS0cos((j)-(t)^)+^^cos0sin(({)-4)^) 

-i:^cos0„sin0  +L^cos(^-^J  -Z^sin0„sin((f)  -^^)] 

+4>[^4,cos((j)  -(})^)  -Z,sin0gSin((|)  -4)  J  (55) 

-ZjSin0gCos0cos((i)  -(j)^)  +Z,cos0gSin0 


The  Fourier  expansion  of  eq.  (55)  is: 


ikr 

{0  [i/j(A:p^sin0)(cos0sin0^^^^+Z,^„) 


^  /2COS0  77sin0 

+7  (Arp  sin0)(-^^°^^  K.  -cos0  sin0Z,  - - ^L,  )] 

n\  r  o  /V  T  -  /-I  (p,«  g  t,n  j _ 


A:p  sin0 


Arp  sin0 


+4)  [/y„'(Arp^sin0)(i:^  „  -  cos0sin0/^  J 

+y  (Arp  sin0)(cos0^sin0Z  - — ^ ,7  )]} 
”  "  ^  ytp^sin0  ytp^sin0 


(56) 


where  J^(x)  is  the  Bessel  function  of  the  first  kind  of  order  n,  (x)  is  the  derivative  of  j  (x) 

with  respect  to  x. 

For  the  scattered  field  in  the  0  =  0  direction,  eq.  (55)  is  further  simplified  and 


transformed  to: 


(57) 
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m.  ZERO  ON-AXIS  BACKSCATTERING  CROSS  SECTION 


A.  INTEGRODEFFERENTIAL  EQUATIONS  FOR  THE  SURFACE  CURRENTS 
The  impedance  boundary  condition  links  the  tangential  components  of  the  total 
electric  field  £  to  the  tangential  components  of  the  total  magnetic  field  ^  on  a  surface 
through  its  surface  impedance,  which  is  a  function  of  the  electromagnetic  properties  of  the 
material  of  the  surface.  This  condition  is: 

X  X  =  2"  («-  X  H^)  (60) 

where  ^  E'”‘^  ^nd  are  the  total  e  ^nd  ff  fields  and  Z and 

Z  ■  are  the  normalized  surface  impedances  on  5  +  and  S  ^  respectively.  Since  only  the 
tangential  components  of  and  are  involved  in  eq.  (60),  we  shall  use 

tan  Ian 

^mc±^  denote  the  tangential  components  of  the  scattered  and  incident  fields  on  the 

tan  tan 

surfaces  ,5"  ^  and  S  ^  ■ 

In  terms  of  surface  electric  and  magnetic  current  distributions  eq.  (60)  becomes: 

±  fi^  ^  (61) 

For  anisotropic  materials,  it  is  more  convenient  to  utilize  matrix  notations.  In  eq.  (61)  the 
vectors  and  can  be  considered  as  two-element  column  vectors,  with  each  of  the  (])- 
component  designated  as  element  1  and  the  t-component  as  element  2.  Z"  and  Z '  are  then 
two-by-two  matrices. 

The  cross  product  of  the  unit  outward  normal  n  *  with  such  a  two-element  column 
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vector,  which  is  composed  of  ({)-  and  t  -components,  can  be  represented  by  the  following 
matrix  multiplication: 


L  =  V  V 


(62) 


where  F 


0  -1 

1  0 


F 


T  _ 


-F  and  F'^  =  -/.  Hence,  in  terms  of  the  matrix  F,  the 


impedance  boundary  condition,  eq.  (61),  can  be  written  as: 


Z-K~  =  -FL 


(63) 


It  should  be  noted  that  ^  and  x  fi=  involves  and  which 

can  be  computed  from  ^  +  L~  viaeqs.  (51)  and  (52).  Since  the 

difference  currents  g*  -  g~  and  g^  _  g~  can  also  be  obtained  from  eqs.  (51)  and  (52), 
the  impedance  boundary  condition,  eq.  (63),  if  written  as  a  relation  between  the  sum  and 
difference  currents,  becomes  a  set  of  integrodifferential  equations  for  the  sum  currents.  Since 


Z^  K*  =  1  Z"  [  ^  +  (K*-  K~)  ] 


(64) 


Z-  K  =  Iz-  [K-{K-K)] 


(65) 


eq.  (63)  can  be  written  as: 


(67) 


V  r  =  -  1  Z-  [  f  -  {K^-  K')  ] 

Define  Z  and  A  as; 

Z  =  1  (Z-  t  Z-)  (6S) 

A  =  i  (Z-  -  Z  )  (69) 

By  subtracting  eq.  (67)  from  eq.  (66)  and  substituting  Z  and  A  in  the  result,  we  have 

Z  (i "  -  Z  ■)  =  Z  £  +  A  -  ^■)  (70) 

By  adding  the  two  equations  together,  we  have 

V  L  =  Z  {K^  -  K')  +  A  K  (71) 

Assume  that  Z  exists,  we  obtain: 

K*  -  K~  =  Z-^V  L  -  Z-^AK  (72) 

Substituting  eq.  (72)  this  into  eq.  (71),  we  have: 

V  -  L  ^  -V  [Z  -  AZ-^A]K  -  V  AZ-^VL  (73) 

Since 

±ZZ^  =  X  (r  X  „-)  =  =  (£^  (74) 
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and 


±VK-  =  =  -H 


Eqs.  (72)  and  (73)  can  be  written,  with  the  understanding  that  and 


tan  tan 


H  '  +77  =2/7  ’  as: 

tan  ■‘tan  '^'‘tan 


C  -  E-  +  2EZ  =  [Z  -  AZ-^A]/:  +  ^Z-^VL 


-  VZ-^A  K  -  VZ-^VL 


As  r  approaches  froni  outside  S,  the  tangential  components  of  £  and^ 
are  derived  in  Appendix  C  as  given  by  eqs.  (51)  and  (52).  They  are; 


Er(M=^-E^ 


f^/^y^^’''*APo(0[^c|,cos((t)-(^^)-/:^sin((j)-c|)Jsine„(5„)]G(r-rJ 


i  d  r  I  r^Tt 


- - ^  f  f  as  d6[X^-^G(f-r  )  +  p(s)K  —G(r-r  )] 

4Tcp^(5)a(l)Joio  ''oil 

■  sin(4)-(l)^)cos0/5) 

+/,cos((()  -(i)^)cose^(5)sin0^(5^)  -/^sin0g(.s)cos0„(5^)]G(r -f^) 

+/,COS0  (5)cos0  (5j  +/*sin((i)  -(1) Jsin0  r5)]G(f-f  1 
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*sc± 


E,  ((|),5)  =  t  ■  E 

=  ^  /^'|“*Ap,(s„)[K*sin{(l)-<t.„)sme/s) 

+ Kpos{(^-<^Jsmd^(s)sind^(sJ  +K^cosd^(s)cosd(sJ]G{r-fJ 


471  55^0  Jo 

k  d  c  ^  c'^T^ 


ad)  ^  ^a^ 

~  o  o 

Lj  sme^(5)  0080^(5^)  +Z^cos(4)  -(j)^)  0080^(5)  Bind  ^{sJ]G(r-rJ 
COS0^(5)  d 


(79) 


k 

4% 


Pg(5)  dn 


15* 


[^4.  C08((|)-(l)^) 


where 


dn 


-  L^sin((p-<l)Jsmd^(sJ]G(r-rJ 

is  the  normal  derivative  taken  in  the  limit  as  r  approaches  the  surface  S"'  or  S 


respectively  from  outside  of  S.  It  is  evaluated  as  a  limiting  value  in  this  report  and  should  not 
be  confused  with  the  Fourier  index  n.  Similarly,  we  have: 

^  |^y^^"^/(l)^p/5^)[Z^cos((j)-4)J-i:^sin((l)-(j)^)sin0^(5^)]G(f-f  ) 

i  d  r  ^  r2n 


47ip^(s)a(t)JoJo  "  *'a4>^  ' 


^  sm(<i>-<|>.)cos9/5) 

+^,cos((t)-(})^)cos0  (5)sin0  (5j-i:/(})^,5^)sin0  (5)cos0  (5^)]G(f-fJ 


(80) 


4n: 


k  (5)- — 
'  ^  dn 


[  f  (5^)[i:cos(({)-4)  )sin0„(5)sin0^(5  ) 

.±  ^0  Jo  ^  ^ 


+K^cosd  (s)cosd  (sj  +E^sm((p -(pJsinO  (s)]G(r-FJ 
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+Z,cos((|)  -(|)^)sin0g(5)sineg(5^)  +Z,cos0„(:y)cose(5^)]G(f-f J 

[Z ,  —G(r-f  )  +  p(s)L  —G(f-r  )] 
47:05^0  Jo  °  °  ‘‘’06  °  ^  '0.V  ^ 


471  dsJo  Jo 

k  d  r  ^  r2n 


s  ^ 

-  K,sme^{s)cosQ^(sJ  +Kfos(<i>-<^Jcosep)sme^(sJ]G(r-rJ 


47t 


COS0  r^)  3 


p  r^)  dn 


///_,’**  Ap  A)  cos(<)>-4>J 


-  i:^in((i)-4)^)sin0  (5j]G(f-rJ 


(81) 


As  r  approaches  S*  and  S  ,  p-’p^(5)  and  z-'zj^s),  the  Fourier  expansions  of  eqs.  (78)  and 
(79)  are: 


.P 


PXs) 


1  rl 


4  05 


ik 
4 

k  d  rl 


K„(^)  - 


dn 


/y^oPgi^o)  sin0g(5)  (G„^j-G„.j)  „ 


2  ^ [smd/s)cosd/sJ G„ 


-  -  COS0/5)  sin0/5^)(G„,j  +  G„-i)] 


K  (5)  - - 

dn 


.  /J^oPg(‘^o)[cOS0^(5)  COS0g(5^)  G„ 


1  . 


-  sin0/5)  sin0/5^)  (G„^,  +  G„.,)]Z^^„ 


(82) 
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e:t(^  = 


{  is  />^  PA)  [A  Gj 

ik^  rl 


f  ^o  Pg(V  [(^osd^(s)cosd^(sj  G„ 

J  O 


1  . 


+  -sine^(5)sin0/5j(G„,i  +  G^_,)] 

cosQXs)  d 


p  ^(5)  dn 


f‘ds^  P/O  "  G«-i)  h.n 

d  O 


k 
4 

kncosQis)  ., 

- ^  />»  P/O  (G..,  -  G..,)  i,  . 


ik 


4pg(^) 

cos  0(5)  s 


p  G)  dn 


ikn  r  i 


fds,  p/5j  sin0  (5^)(G„,j  -  G„_i) 

d  o 


2pXs) 


Pg(0  i\  cos0^(5)sin0/5^)(G„,j+G„_^) 


■sin0^(5)cos0^(5^)  GJ 


(83) 


Eqs.  (80)  and  (81)  lead  to  similar  forms  for  and  by  replacing  £  with  _  jf 
and  £  with  £  in  eqs.  (82)  and  (83).  Therefore, 


O)  *  c  w "  Af,  4  (5) 


(84) 


where  the  elements  of  the  two-by-two  matrices  and  are  integrodifferential  operators 

on  the  elements  of  the  column  vectors  £  ,  £  defined  by; 
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n  d 


KMs)  =-/>AW[f  “"9/0(g.-.-g,..)+-^^gj  kjo 


PjW 


M 


2 

^  C  ^  /^  TJ'  /  \  k  /'  / 


M 


naih,rk^)  = 

,.2.A,W  = '  £  /  'A  p/A  Gj 


•'  o 

2 


(85) 


P/O  [cos0^(5)cose(5jG„ 

o 


1  . 


+  -sine^(5)sine^(5^)(G„,i  +  G^_,)]  KJsJ 


K,nhJ^)  =  -f-^f  '^.PA)cos0g(5)(G„,j  -G„.j)Z^„(^J 


i/c 

T 


2  05. 

2k  (5)-^ 
5  r  / 


0^2 


f^^oP/O  sin0/^)  (G„,j  -G„_i)Z^  „(0 


/^  ^oPg(0  [sin  QgC-y)  cos0g(^^)  G„ 

-  j  cosd/s)  sind/sj  (G„,  j  ^G„J]LJsJ 


2k  (5)-  — 
^  d/7 

1  . 


dn 


/J^oP/O  [  COS0^(5)  COS0^(5^)  G„ 


+  -  smd/s)sme^(sj  (G„,i  +G„.i)]Z^^„(5^) 


(86) 


k 

2cos0^(5)  ^  5 

0 

4-  _ 

r 

4 

Pg(^)  0« 

s-  <3" 

S"_ 

Jo 

T  -  ^  O 


A77COS0  (5)  ^  , 

- ^f  dso(s)(G  ,-G  ,)^a(5) 

2p  (5)  Jo  ^  (t),«K  O^ 


ik 


2cos0g(5)  ^  a 

0 

+ 

r 

Pg(-5)  0« 

5.  dn 

f  '^■^<,Pg(Osin0/O(^«.i  -^„-i)4,«(0 

T-  d  O 


ikn  r  1 


/^^.PA)[2 '''’s0g(^)sineg(O(G„.i  +G„.i) 


-  sin0  (5)  COS0  J5^)  G J  L^  „(5^) 


Similarly,  we  can  write: 
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(87) 


Hr(s)  +  Hn^)  =  -  N  K  (s)  -  M  L  (5) 

tan,??''  ^  n  n  ^  ^  n  n  ^  ^ 


Combining  these  with  Fourier  components  of  eqs.  (76)  and  (77),  we  obtain  the 
equation  for  the  sum  currents  as: 


< 

K 

ri- 

Z-AZ^^A 

AZ-^V 

=  2 

inc 

tan,?? 

N 

n 

M 

Yl 

<1 

-vz-^v 

L 

L 

n 

Uji,n 

(88) 


Note  that  given  the  incident  wave  on  the  shell,  this  is  a  system  of  linear  integrodifferential 
equations  of  the  second  kind  for  the  unknown  currents  if  ,  T  .  It  is  assumed  that  unique 

n  n 

solutions  to  eq.  (88)  exist  for  all  n. 


B.  ZERO  ON-AXIS  BACKSCATTEREVG 


Since  the  shell  is  rotationally  symmetric,  a  wave  incident  along  the  axis  of  the  shell 
(the  z-axis)  can  be  assumed  to  be  linearly  polarized  in  the  y-direction  without  loss  of 
generality  (see  Figure  2).  Let  the  fields  on  the  surface  of  the  shell  be: 


E  (({>,5)  =  j;  e  « 

-Frinc.t  .  ^  -ikz^ 

H  ((}),5)  =  X  e  * 


(89) 


which  represents  a  plane  wave  propagating  in  the  -z  direction.  Because 
y  =  cos4)4>  +  sin4) p  and  p  =  sine^f  +  cose^ri+ .  then  they  become 
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£'"‘'(4),  5)  =  [  coscfxf)  +  sincj)  (sin0^^  +  cosO^t?")  ]  e 

o 

i7*"‘^(4),  5)  =  [-sin<|)4)  +  coscj)  (sin0  t  +  cos0„«"')  ]  e 


The  tangential  components  for  the  incident  fields  on  the  surface  of  the  shell  are  expressed  as: 


7--  me 

e;"‘ 

-rrinc 

^4> 

Rinc 

sincf)  sin0^ 

-sine}) 
coscj)  sin0^ 


-ikz 

e  ^ 


The  Fourier  coefficients  are  nonzero  only  for  n  =  ±  1 : 


1! 

1  -ikz^ 

~e  ^ 

2 

7—'  me 

7-  me 

sin0  e  = 

2  - 

inc 

-  ^r.-l 

= 

1  -ikz^ 

2 

TJ 

h:t  = 

— sin0  e  - 

2  ^ 

Kl 

Because  only  the  n  -  ±\  terms  in  the  incident  field  are  nonzero,  only  and 
can  exist.  Therefore,  only  the  77  =  ±1  components  of  eq.  (88)  have  to  be  considered: 


Mbi  -^±1  Z-l^Z-^b.  AZ-*F  I  Kx 
^±x  HiJ  [  vz-^b.  -vz-^vj  J  4 


The  incident  fields  in  eq.  (93)  can  be  written  in  the  form 
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^  me 
'tan,±l 

inc 


2 


1 

F  /  sin0^ 

* 

±  i 
sin0„ 


e  * 


(94) 


which  satisfies  the  following  equation; 


0 

-il 

^  me 

=  ± 

^tan,±l 

il 

0 

0inc 

^tan,±l 

^mc 

^tan,±l 

(95) 


or,  equivalently, 


0  -il 
il  0 


H. 


?  me 
taruil 

ine 


taD,±l 


(96) 


Therefore,  if  both  sides  of  eq.  (93)  are  multiplied  with  the  factor  j  j 
right  hand  side  of  eq.  (93)  is  equal  to  zero  due  to  eq.  (96): 


0  -il 
il  0 


Hi 

0  -il 
il  0 


Z-AZ'^A  AZ'^F 

VZ'^^  -vz-^v\ 


K 


±1 


L 


±1 


1, 

1 

O 

I 

^tan,±l 

1  + 

[ 

il  0 

J 

r}inc 

then  the 


(97) 


=  0 


From  eqs.  (85)  and  (86) 
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//IK  K  -// 


M.J  ^A/„  0 


Therefore,  if  the  matrix 
becomes: 


Z-A  Z-^A  A  Z'^Z 

also  commutes  with 

0  -il 

ZZ'^A  -ZZ'^Z 

il  0 

,  eq,  (97) 


Ki  -^±i]  Z-AZ“^A  AZ“'Z 


ZZ-'A  -ZZ'^Z 


0  -il]  1 

1  +  (  =0  (99) 

k  oJ  r„ 


The  existence  and  uniqueness  of  solutions  to  eq.  (97)  implies  that 


0  -il  ^.1 
1  +  (  =  0 
il  oJJ  A., 


(100) 


or,  equivalently 


±  iZ/^j  =  0 


(101) 


Substituting  this  result  into  eq.  (59),  we  conclude  that  =  0  if  the  matrix 


Z-AZ'^A  AZ'V  .  [o  -z/1 

commutes  with  .  This  is  a  suflBcient  condition  to  eliminate  on- 

VZ-^A  -V  Z'^V  [il  0 


axis  backscattering. 


IV.  IMPEDANCE  MATRICES  FOR  ZERO  ON-AXIS  BACKSCATTERING 


A.  ELEMENTS  OF  Z  AND  A 


Z-AZ'^A 

AZ-*F 

0  -il 

It  can  be  verified  that  the  matrix 

commutes  with 

FZ‘‘A 

-V  Z~^V 

il  0 

and  only  if; 


VZ-^A  =  -AZ-^V  (102) 

Z+FZ-‘F  =AZ-iA  (103) 


where  both  Z  and  A  are  2  x  2  matrices.  Under  the  assumption  that  the  inverse  of  Z  exists, 
we  analyze  eqs.  (102)  and  (103)  as  follows: 

Because  the  determinant  of  the  two  sides  of  the  eq.  (102)  has  opposite  signs,  we  get 
detZ'MetA  =  0  - Since  Z  is  assumed  to  be  invertible,  det  A  =  0. 

Because  of  the  identity: 


z-i  =  -  -J—VZ^V 

detZ 


(104) 


and,  by  multiplying  V  to  both  sides  of  eq.  (102); 

Z-^A  =  V  A  Z-^V  (105) 


eq.  (103)  can  be  transformed  to 


A(VAZ~'^V)-V 
-AVAV(VZ-^V)  -  VZ-^V 

1 


detZ 


[l+(AVy]Z- 


where 

^12(^12  ~^2l)  ® 

0  Ai(^12-^: 

In  the  last  equation,  we  utilized  the  fact  that  AjjA22  =  Aj2A2j  because 
Substituting  eq.  (107)  into  eq.  (106)  we  conclude  that; 

1.  If  Zjj  #  0  or  Z22  *  0,  then 

1  ^^12(^12-^21)  = 

1  -  A2j(Aj2-A2j)  =  detZ 
Z  =  Z^ 

2.  IfZj,  =  Z22  =  0,  then  detZ  =  -Z^^Z^^  *  0 

1  +  Aj2(Aj2~A2j)  =  “  ■Z’22 

1  -  A22(Aj2~A2j)  =  ~  Z2J 

On  the  other  hand,  substituting  eq.  (104)  into  eq.  (102)  yields: 


(AF  y=AFAV= 


^11^22  ^12 


0 


^11^22  ^^21 


(106) 

(107) 

1) 

detA  =  0 . 

(108) 

(109) 

(110) 

(111) 

(112) 
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(-^12  ^21)^11 


(113) 


•^22(^12  ^21)  (^12  ^21)^22 


(114) 


■^11^22  "''  ■^22'^11  ^■^21^12  ^■^12'^21 


(115) 


For  Case  1,  and  ^re  not  both  zero  and  Z  is  symmetric.  Since  Z^2  =  Z2j ,  eq. 
(1 13)  or  (1 14)  requires  that  Aj2  =  A2J .  Therefore  A  is  also  symmetric.  Eqs.  (108)  and  (109) 
both  requires  that  detZ  =  1 .  Multiply  eq.  (115)  with  either  A^^  or  A22  and  replace  AjjA22 
with  A  2  we  find  either 

iA,2’ 

■^22^11  ~  ^^12^1l'^12  ■^11^12  ~  ^  (11^) 


or 


^11^22  ^•^12^22'^12  -^22^12  ~  ^ 


(117) 


Eq.  (1 16)  can  be  solved  for  A^,  in  terms  of  Aj2 ; 


^  _  ■^12^V‘^12  ■^11-^22  ^ 

'^12  “ 


■^12  ^  ^ 


‘11 


Z 


‘12 


'22 


^22 


Multiplying  to  both  sides  and  replacing  Aj^A22  with  ^2^  again  result  in: 


A  - 


Z,-  ±  i 


"22 
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Hence  both  Aj^  and  can  be  given  in  terms  of  A22  as: 


Similarly,  from  eq.  (117); 


A„  = 


/  7  ±,•^2 
^12  ' 


Ai,  = 


7 

V  ^22  j 

(  Zj2  ±  i 


A 


22 


Z 


A 


22 


22  J 


^22 


= 


11  J 


(118) 


(119) 


Eqs.  (118)  and  (1 19)  are  equivalent  and  both  include  the  situation  A  =  0 .  The  particular  one 
to  use  is  a  matter  of  convenience,  especially  when  either  Z^j  =  0  or  Z22  =  0.  Note  that 
Zj2  =  Z2J  =  ±i  if  ZjjZ22  =  0  because  detZ  =  1. 

ForCase2,  Z^j  =  Z22  =  0  and  Zj2Z2j  O.Eq.  (115)  requires  that  A  ^2  =  =  0, 

therefore  Ajj  A22  =  Aj2  A2J  =  0.  Eqs.  (1 11)  and  (1 12)  reduce  to: 

^12  =  =  -  1  (120) 


^  ^12  ~  ^21  ~  ^  no  further  requirement  on  A  except  that  at  least  one  of  A^  j 

or  A22  vanishes.  If  Zj2  =  -Z21  =  ±i,  then  eqs.  (113)  and  (114)  requires  that 
A]i  =  A22  =  0.  Hence 
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=  Z-  -  2  = 


(121) 


0  ±i 
+i  0 

This  is  the  only  situation  when  detZ  *  1 . 

B.  Z'ANDZ 

The  impedance  matrices  for  zero  on-axis  backscattering  specified  in  Sec.  A  above  can 
be  stated  in  terms  of  Z  ^  and  Z  ‘  directly.  When  z "  and  Z  “  are  skew-symmetric,  the  only 
possibility  is  eq.  (121).  When  Z*  and  Z"  are  symmetric,  the  requirements  are: 


detZ""  =  1 

(122) 

detZ'  =  1 

(123) 

det(Z"  -  Z')  =  0 

(124) 

Eq.  (124)  is  just  detA  -  0.  To  show  that  conditions  (122)  through  (124)  are 
equivalent  to  eqs.  (108)  to  (115),  note  that  eqs.  (1 10),  (1 13),  and  (114)  are  satisfied  trivially. 
Eqs.  (108),  (109),  (111)  and  (112)  become: 

detZ  =  1  (125) 

and  eq.  (115)  becomes: 
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(126) 


■^11^22  ■^22^11  ^'^12^12 


Since  Z"  ^  Z  +  Z~  =  Z  -  L,  eqs.  (122)  and  (123)  are: 

(^11  +^ll)  (^22  +  ^22)  "  (^12  ■^^12)^  =  ^  (127) 

(Zjj  “Ajj)  (Z22  “  A22)  ~  (Zj2  “  Aj2)^  =  1  (128) 

Except  for  an  overall  factor  of  2,  the  difference  of  eqs.  (127)  and  (128)  is  eq,  (126);  the  sum 
of  these  two  is: 

■^11  ■^22  '''  '^ll  “^22  ~  -^12  “  Aj2  =  1  (129) 


which,  when  combined  with  detA  =  0,  is  just  detZ  =  1 .  Together  with  eq.  (126), 


det(Z  +  aA)  =  detZ  +  adet 


^12 
'■^12  ^22 


+  cc  det 


^11  ■^12 
^12 


+  cc^detA 


=  detZ  =  1 

for  any  a  Therefore,  with  a  =  ±  1 ,  detZ "  -  detZ  “  =  1  •  Hence,  with  symmetric  Z "  and 
Z  Z  conditions  (122)  to  (124)  specify  completely  the  impedance  matrices  which  eliminate  on- 
axis  backscattering. 

It  should  be  noted  that,  if  the  shell  is  a  closed  surface,  then  the  impedance  boundary 
condition  closes  off  the  inside  of  the  shell  from  outside.  Therefore,  only  detZ'^  =  1  is 
required  to  eliminate  on-axis  backscattering.  This  is  an  extension  of  Weston’s  work  [7]  to 
anisotropic  impedance  coated  body  of  revolution. 
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V.  CONCLUSIONS 


A  sufficient  condition  for  eliminating  the  on-axis  backscattering  of  an  anisotropic 
impedance  coated  shell  of  revolution  has  been  deduced.  The  outside  and  inside  normalized 
surface  impedances  Z  ^  and  Z '  with  which  this  sufficient  condition  can  be  satisfied  have  all 
been  found.  One  exceptional  situation  is  when  the  impedance  matrices  are  equal  and  skew- 
symmetric  with  their  determinants  equal  to  -  1 .  All  other  cases  require  that  the  two  matrices 
be  symmetric,  their  determinants  be  unity,  and  the  determinant  of  their  difference  be  zero. 

Based  on  the  duality  of  Maxwell  equations,  Weston  [7]  argued  that  an  object  coated 
with  a  unit  normalized  surface  impedance  will  have  zero  on-axis  backscattering  if  the  object 
is  invariant  under  a  90°  rotation  around  its  axis  of  symmetry.  It  is  clear  from  Chapter  III  of 
this  report  that  the  duality  property  is  crucial  in  determining  the  sufficient  condition  for  zero 
on-axis  backscattering.  Results  of  this  work  therefore  should  not  be  limited  to  shells  of 
revolution  only.  This  aspect  of  the  problem  is  being  completed  and  will  be  published  shortly. 

The  shell  under  consideration  can  be  a  closed  one.  A  body  of  revolution  coated  on  the 
outside  with  an  anisotropic  surface  impedance  Z  ^  can  be  modeled  as  a  shell  having  Z^  =  Z- 
Since  the  impedance  boundary  condition  separates  completely  the  exterior  of  the  body  from 
its  interior,  only  the  conditions  on  Z*  need  to  apply,  i.e.,  Z"  rnust  be  either  symmetric  or 
skew-symmetric,  with  detZ^  =  ±1-  This  is  an  extension  of  Weston's  result  to  anisotropically 
coated  bodies. 

It  should  be  noted  that  a  properly  shaped  and  coated  body  can  always  be  carefiilly 
deformed  into  a  shell  without  generating  any  on-axis  backscattering.  This  procedure  yields 
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the  condition  Z  ^  =  Z  ‘  By  beginning  with  a  shell,  we  are  able  to  relax  this  condition  to 
det(z "  -  Z " )  =  0  for  the  symmetric  case. 

Results  of  this  work  make  available  a  wide  class  of  models  which  must  have  zero  on- 
axis  backscattering  cross  section.  All  general  purpose  numerical  codes  for  computing  the 
scattering  cross  sections  of  anisotropic  impedance  coated  objects  should  be  checked  for  their 
accuracy  against  a  selected  group  of  such  models.  Such  comparisons  should  provide 
indications  of  an  error  bound  of  the  particular  algorithm. 
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APPENDIX  A:  VECTOR  CALCULUS  ON  A  SURFACE  OF  REVOLUTION 


Vector  calculus  near  a  surface  of  revolution  will  be  presented  in  this  appendix.  In 
particular,  the  operations  of  the  divergence  and  curl  on  a  vector  field  will  be  given  explicitly. 
The  derivation  will  be  based  on  vector  calculus  in  the  familiar  cylindrical  coordinates 
(P, 

The  axis  of  revolution  is  taken  to  be  the  z-  axis.  In  the  half  plane  of  a  fixed  <|) ,  the 
generating  curve  of  the  surface  is  parameterized  in  its  arc  length  5  with  the  functions  z^(5')  and 

Ks)  =  2g(s)z  +  p^(5)p  (130) 

The  surface  is  thus  parameterized  in  terms  of  (4*,  ‘^)  •  This  generating  curve  may  be  closed, 
open,  discontinuous  or  self-intercept.  Nevertheless,  at  every  point  in  a  smooth  segment  of  the 
generating  curve,  an  outward  unit  normal  to  the  “outer”  surface,  denoted  by  /r  can  always 
be  defined.  The  parameter  5  can  be  oriented  so  that  the  tangent  vector  i  to  the  generating 
curve  satisfies  the  relation: 

«"  =  4>  X  f  (131) 


Note  that: 


—  -  COS0  (5)f  +  sine„(5)p 
as  ^ 


(132) 


where  cos0^(j’) 


sin0g(v) 


—  Po(‘^)i  is  the  angle  measured  from  z  to  i , 
ds  ° 
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defined  to  be  a  continuous  function  of  5  along  any  smooth  segment  of  the  generating  curve. 
Since 


d^f 


d 


=  (-sin0  (5)z  +  cose,(5)p)— 0^(5)  =  K„(5)/7 
ds^  °  da  ^ 


(133) 


where  k  (5)  -  —  0„(‘5')  is  the  signed  curvature  [8]  of  the  generating  curve:  k  is  positive 
when  0^  is  increasing  with  s  and  1 1  /  |  is  the  radius  of  curvature  of  the  generating  curve. 

In  the  neighborhood  of  the  surface  of  revolution,  the  space  can  be  described  with  the 
orthogonal  coordinates  (cj),  5, «),  n  -  0  being  the  surface  and  n  >  0  going  in  the 
direction.  In  relation  to  the  cylindrical  coordinate  system,  ({)  is  identical  and,  for  «  0,  the 

transformation  of  dz,  dp  into  ds,  dn  is 


(134) 


ds 

COS0 

s 

sin0„ 

dz 

dn 

-  sin0 

a 

COS0 

sj 

dp 

The  transformation  of  ds,  dn  into  dz,  dp  follows  the  same  rule  as  the  matrix  is  unitary.  This 
equation  also  gives,  for  «  0 ; 

05'  dn  _  dz  _  dp  _ 

dz  0p  ds  dn  ° 

^  _  0^  _  0P  _  _  0-2^  _  sin0 

0p  0z  05  dn  ° 


For  a  vector  field  A  -  ^  +  A^  t  +  A^  n^ ,  A^  and  A^  transforms  fi'om  A^  and  A^ 

according  to  eq.  (126).  On  the  other  hand. 
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(136) 


8 

=  ^ 

8 

COS0 

8 

-  sin0 

8 

8z 

8z 

85 

dz 

dn 

g 

85 

g 

dn 

8 

85 

8 

dn 

0 

-  sin0„ 

8 

+  COS0 

8 

8p  ■ 

8p 

85 

dn 

g 

85 

g 

dn 

Substitute  the  above  into  the  expressions  in  the  cylindrical  coordinate  system  result  in  the 
following  as  n  ^  0 : 


V  ■  A  = 


p  84)  ** 


=  -[ 


P,  d4>  ^ 
P,  54) 


d  .  1  8  ,  ,  s 

— + - (p^j 

8z  ^  p8p  P 

ds  dn° 


sin0„ 
( - ^ 


8  cos6^ 

—)A,  +  ( - ^ 

5^  '  P. 


dn 


K 


(137) 


(V  X  A).  =  —A^  -  —A  = 
8z  P  8p  " 


8  .  .  8 .  - 
— A  +  (k  -  — )A, 

85  "  °  8«  ' 


(138) 


cos0„  3  1  ^ 


dz 


(139) 


cos0 

=  ( - ^ 


,  8  .  ,  18, 

+  — )A,  — - A 

dn  ^  P„  84)  " 
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(V  X  A)^ 


sin0„  p  ^  1  ^ 

-  - — (P^*)  “  — ^o]  +  COS0  [ - A 

p  34)  P-*  -  pacj)  ^ 

f  [^-4,  -  |-(P-4t)] 

P.  5(1)  ds  'P 


dz 


■4J 


±±A,  - 

Pg  54)  p„ 


ds  “* 


(140) 


As  a  check,  the  substitutions  A=-Af.,A  =  A  ,  p  =  rsin0,  6  =  9  -  — 

r  f7  rg  ^  g  2  ' 

1  d  _  i  d  d  d 

,  - - -  — ,  —  =  —  result  m  the  correct  expressions  for  V  ■  ^4  and 

r  ds  r  dd  dn  dr 

V  X  ^  in  the  spherical  coordinate  system  (r,  0,  ({))- 
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APPENDIX  B:  END-CAP  CONTRIBUTIONS  TO  THE  SCATTERED  FIELDS 


In  eqs.  (45)  and  (46),  the  surface  S,  consisting  of  the  four  non-overlapping  surfaces 
Si,  S^,  S-i  and  .S'2(Fig.  3),  can  be  shrunk  toward  the  shell.  So  that  S  may  be  the  surface  of 
constant  distance  6  from  the  shell.  In  this  appendix  the  fact  that  contributions  from  the  top 
cap  Si  and  the  bottom  cap  S^  to  the  scattered  fields  are  zero  in  the  limit  that  the  thickness  of 
the  shell  is  infinitesimal  will  be  demonstrated. 

An  example  of  the  geometry  of  the  end  cap  is  shown  in  Fig.  4.  Near  the  end  at  s  =  I, 

define: 


45 


the  normal  direction  at  the  end  as: 

the  tangent  direction  at  the  end  as 

and  the  tilt  angle  at  the  end  as;  a  =  lim  0^  (5),  which  is  the  angle  measured  from  z- 

s~r 

axis  to  The  cross  section  of  the  end-cap  at  ^  =  /  in  a  half  plane  of  fixed  azimuth  angle  (j) 
is  a  semi-circle  with  a  radius  6 .  Note  that  6,  a  constant  distance  distance  between  S  and  the 
shell,  is  a  fixed,  positive  number.  In  Fig.  4,  the  angle  (j)^  is  measured  from  the  axis  n^.  In 
terms  of  the  angles  (j)^  and  a,  the  end  cap  Si  can  be  parameterized  in  the  cylindrical 

coordinates  as ( p,  (j),  z),  where  for  o<(j)^<7r, 

P  =  Pg(0  +  ^  cos(<j)'  -  a) 
z  =  Zg(/)  +  6  sin(<J)^  -  a) 
and  0  <  (j)< 2tc . 

Note  that  on  the  cap  surface,  any  field  vector  may  be  decompossed  to  the  three  directional 
components:  (f)-  component,  (|)^-component  and  n-component,  where  n  denotes  the  outward 
normal  to  5'^. 

In  eq.  (45)  and  (46),  the  contributions  to  the  scattered  fields  from  surface  integrals 
over  the  end  cap^”;  are: 

/j  =  r  [fix  E(r^)]G(f-fJda^ 

JS, 

(141) 

=  Pp"'  [fi  X  E(rJ]G(r-fJ  [pj/)  +  6cos((})^-a)]  6  d(^d(^' 

•I  0  0 


fi,  =  lim  fi  (s) , 

s-*r 


A  A 

ti  =  lim  t(5), 

s^r 
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(142) 


A,  =  r  [« X  )]G(f-f ) 

•'Si 


da 


0/"  [/?  X  H{rJ\G{r-r^)  [  p//)  +  6cos((t)'-a)]  6 


^1=1  [n  ^  da^ 


=  0o'"  [p^(/)  +  6cos((i)^-a)]  6  ^^(})/ct)' 


(143) 


=  f  [nx  E(01*^oG(^<)  dUo 

Js, 

=  f'f"'  {n  X  %)]*V„G(r-f„)  [p,(/)  -X  6cos((t)'-a)]  6  rf,t>/<t>' 

JO  -^0  ^ 


(144) 


From  the  edge  conditions  [5],  on  the  end  cap  S^, 

nxE  =  E^^' where  |£^|  ^  and  <  £26'*^^ 

and  where  and  \H^,\  <  , 

where  E^,  E^,  H^,  are  positive  constants.  Therefore, 


nxE\  <  +£2^ 

nxH\  <  6 -'0^62  +hI 


(145) 


I 


Note  that  f  lies  away  from  the  shell.  If  the  distance  between  f  and  the  shell  is  6  «  —  6  can 

^  k'  ^ 

be  choose  such  that  6  <  ~  6^,  so  that  f  always  lies  outside  with  . 

2  2 

I  I  2  1  C 

Therefore  \G{f~f^\  <  — —  -  and  -7|V^G(f-f)|  <  -  =  Q  for  some  positive 

ko  .  0  ^ 


2  ^2 


k^b 
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constants  C,  Cj  and 

Substitute  these  into  eq.  (141)  and  (142),  and  write  them  as: 


M"*!  I  S  [pXI)  +  6  cos((|)^  -  a)]  [n  x  %)]  G(f-f ) 

|pg(0+Scos((})^-a)|  |£^4)'-£^4)|  |G(r-f)| 
^ ^  ^ ^  6cos((t)'  -  a) I  +E2  )  Cj 

s  2x2  ,  5J  C, 


Substituting  fi  x  fj  into  eq.  (142),  we  have: 


|X 


fj^^d<t>d4>'5  [p//)+5cos((l)^-a)]  [nx#(F)]  G(r-rJ 


<  271^  (p//)  +  6)  ^ Hi]  C, 


(146) 


(147) 


Similarly,  for  eqs.  (143)  and  (144),  we  have: 

|i|;J  <  27r2(p^(/)  +6)  6^  ^E^]  (148) 

IfJ  ^  27t2 (p/O  +  6)  (51/2 (J49) 

therefore  the  limits  of  |^j  ,  |^^|,  ji|;j|,  jij/^l  vanishe  as  6  -  0 . 

Hence,  in  the  limit  that  the  thickness  of  the  shell  is  infinitesimal,  the  contributions  of 
the  integrals  over  the  top  cap  to  the  scattered  fields  are  zero.  Because  the  geometry  of  the 
bottom  cap  is  similar  with  the  top  cap  S^,  the  results  for  are  also  valid  for  S^.  The 
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contribution  from  the  bottom  cap  to  the  scattered  fields  is  also  zero.  We  conclude  that  the 
integrals  in  eqs.  (45)  and  (46)  need  to  be  carried  out  over  S*  and  ^  -  only. 
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APPENDIX  C:  SCATTERED  FIELDS  NEAR  THE  SHELL 


The  object  interested  in  this  report  is  an  arbitrarily  shaped  shell  of  revolution.  Near 

the  surface  of  relvolution,  the  vector  calculus  derived  in  Appendix  A  should  be  used  for  the 

operation  of  the  divergence  and  curl  on  a  vector  field.  In  this  appendix  the  scattered  fields 

near  the  surfaces  S  ^  and  S  “  are  going  to  be  derived  by  using  the  vector  calculus. 

^  dr 

From  Appendix  A,  t  =  —  =  fcos0  (5)  +  p  sin0  (5),  and  =  A  x  ^  The 

ds  ^  ^  ^ 

coordinate  tranformations  between  f  and  f  are' 

O 


=  cos((t)-(j)J 
-f  =  sin(<j)-(j)^)cos0^(5) 
sin((j)-(|)^)sin0g(5) 

=  -sin((j)-4)^)sin0^(5^) 

f  -  f  =  cos0„(5)cos0„(5^)  +  cos((j)-(f)^)sin0„(5)sin0^(5^) 
=  cos(({) -(|)^)cos0„(5)sm0„(5^)  -  sin0„(5) cos0g(5^) 

In  eq.  (51),  the  first  surface  integral  can  be  transformed  as: 


=  -  ^/^yg^’''^4>o^oPg(O{4>[-4sin(4)-<t>o)sin0g(^„)+L^cos((})-(j)J]G 
+  i  [L,cos0^(5)cos0g(5^)  +Z^cos(4)-(|)^)sin0^(5)sin0„(5j 

+L^sin((j) -(t)^)sin0g(5)]  G+n  "[XjCos((j)  -  (j)^)  cos 0„(5)sin  0^(5^) 
-Z,sin0g(5)cos0/5^)  +L^  sin((j)  -  (l)^)cos0^(5)]  G] 


the  second  surface  integral  can  be  transformed  as: 
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=  '*7oX^*‘^'*’”‘*'’‘’**^‘’X'^t^4>  ““<'*’  ■  ‘t’J  -■^,sin(<t>  -<()Jsinaj(iJ]  G 
■K[^j,sin(cl)  -ij>Jsin0^(s)  +A:,  0056^(5)  cos9j(sJ 

cos((l)  -  (f)^)  Sin0,(5)  sin9^(5j]  G+n"[K^  sin({i)  -  ({)^)  0080^(5) 
-r^sin0g(^)cos0^(5^)  +X,cos(cl)  -  cj)  J  cos0^(5)  sin0„(5  J]  G] 

Since  Vij;  =  (4)  — —  —  +  i  —  +  fi*  —)  ij;,  the  third  surface  integral  in  eq.  (51)  becomes- 
p„(5)  04)  05  a« 

o 


X/X'''>apa) 


^  I  r  lTi 

^0  Jo 

'  I  rlTZ 

^0  Jo 


=  -  /„  '!>„  *^,i) 

=  -/v/;rrf4,Ap^(.^ 


=  -z4) 


■lr2n 

>0  Jo 

1  0  rlr2n 


^  0  0  ^ 
— 5 — ^G+K  —  G 

pM  5^0 


S  1  0  ;  0  ..  0  ^ 

(t) - +/  - - 

'’ds^  °0« 


V 


PgC'y)  5<f) 


-G+K^ 


V  Pg(0 

— —G+K^—G^ 

.  p  (5  )  04*  05 

\  ^  O''  o  ) 


ds^ 


— —G+K^—G 

P/O  ds^  j 


G 


oj 


(153) 


From  eqs  (138)  to  (140),  the  curl  of  eq.  (151)  is: 
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"  ~  -I,sin0^(5)cos0^(5j 

+Z,cos((t)-(j)^)cos0g(5)sin0^(5^)]  G 

^  '  Jj )  ~  sin0^(5)  +Z,cos0g(5)  cos0g(5^) 


-kt 


:  i(  cose  (s)  a  ^ 


V  P.(-) 


dn 


+I,cos((l)-(l)^)sin0^(5)sm0^(5^)]  G  } 

/o  Pg(-^o)  [-^4.  cos(4)  -  4)^) 

-Z^sin((})-(|)^)sin0^(5j] 


1  d  /  271 

"p^^/o  /o  -^,sin0^(5)cos0„(^^) 

+Z,cos((f)-4)^)cos0g(5)sin0g(5^)]  G  } 

)1  d  12 

-4>^)sin0^(5)  +Z^cos0^(5)cos0„(5^) 
+Z,  cos((|)  -(|)^)  sin0^(5)  sin0g(5^)]  G 


sin0  (5)  a 
- 1 - + - 

V  Pg(^)  ds 


fo  ■4>o)  ->Z^fSin((|)  -<t)^)sin0g(5  J]G  j 


(154) 


As  r  approaches  from  outside  S,  the  tangential  components  of  as  given  by  eqs. 

(152)  to  (154)  become; 


i 

j 

I 
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'  p.w|-G(-gi 


4tzp  (s)  3(1). 


L^cos((t)-ct)^)cose^(5)sin0^(5^)-I^sin0^(5)cos0g(5^)]G(f-f;) 

I  r2Tz 


(155) 


k 

4tz 


'  ^  dn 


.  fo  fo  ''^°^^oP/^o)t^/Cos((t)-(l)^sin0g(5)sin0^(5^) 

+Z,cos0g(5)cos0g(5^)+Z^sm((t)-(})Jsin0^(5)]G(f-7;) 


£r(4),5)  =  f 

"  ^ /(,Y/’'^o^‘i’oPg(^o)[^<^sin((f)-(|)^)sin0^(5) 


^/:,cos(<{)  -<l)^)sin0g(5)sin0g(5^)  +A:,cos0^(5)cos0(5^)]G(r-f;) 


k  d  r  I  r'^.Ti 

' 

■^fSin0g(5)  cosQ^(sJ  +Z^cos((j)-(j)^)cos0^(5)  sin0g(5^)]G(f-f^) 
|‘cos0g(5)  a 


A: 

47r 


p  (5)  dn 


[^4,  COS((l)-(t)^) 


-  z,sin((|)  -  (|)^)sin0  (5^)]G(f-f  ) 


(156) 


where  — 
dn 


is  the  normal  derivative  taken  in  the  limit  as  f  approaches  the  surface  or  S 


respectively  from  outside  of  5'.  It  is  evaluated  as  a  limiting  value  in  this  report  and  should  not 
be  confused  with  the  Fourier  index  n.  From  the  duality  principle,  the  tangential  components 
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of  can  be  obtained  by  replacing  with  jf  with  and  with 


ik^  r  I 


k  d  r  l  rln 


471  ds 


fofo  sin(<J)-(|)^)cos0^(j) 


-^,cos(<l)-(l)^)cos0  (s)sin0  (5^)-^X<|)^,5^)sin0  (5)cose  (5j]G(r-f;) 


(157) 


47T 


'  ^  dn 


'~  '  g 
'  I  r2n 


fo  ’'^o^<l>oPg(O[^rC0s((|)  -(|)^)sin0g(5)sin0^(5^) 


+/:,cos0g(5)cos0^(5^  +/:^sin((j)  -(|)^)sin0^(5)]G(f-f^) 


ik^  r  l  r2% 


+I-,cos((i)  -(j)^)sin0g(5)sin0^(5^)  +Z,cos0^(5)cos0(5^)]G(f-r^) 

-J-Af^f  ^^dsM  [L.  —G(r-r  )  +  p  (5  )L, —G(r-r  )] 
4T:dsJoJo  ^  ^  ^d(b  ^  ^ 

~o  o 

-  J^,sin0g(5)  cos0g(5^)  +Z,cos((})-4)J  cos0^(5)  sin0g(5j]G(r-fJ 


(158) 


47t 


COS0  (5)  a 


p  r5)  a« 


[K^  cos((j)-(f),) 


-  A:,sin(({)-({)^)sin0  (5^)]G(r-f  ) 


Because  ofthe  rotational  symmetry  of  the  scatterer,  G(r-f^)  =  ‘^"^^niP’^yPo’^J- 

n  =  ~oo 

The  Fourier  expansions  of  eqs.  (156)  and  (157)  are; 
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2Jo  "‘2 
1  r  I  ,  .  .  srk'^ 


P.C-S)  5^0 


f  S  /'*AW  cos6/.)(G...-G..J  L^_, 


ik 

~4 


k(^)  -  — 

«  aw 


sine/5) 


^  I  ^  /J^oP«(0  [sin  0g(^)  cos  0g(5^)  G„ 


Ki^)- 


-cos0/5)sin0^(5j(G„,,  +G„_i)]Z,„ 

r  ^^oP/-^o)[cOS0  (5)COS0  (5jG„ 
s*f 

4^in0g(^)sin0/O(G„.i 


- 1  is  />“  'i- 

o 

+  Pg(0  [cos0^(5)cos0g(5^)G„ 


4  Sin0^(5)sin0/5,)(G„,,.G„., )]/:,,„ 


k 

4 


cosd  (s) 


Pi^) 


dn 


/'*.  P/O  (G,..  +  G,_.)  i^,, 

*/  o 


Ancos0r5)  . , 

-5^  /  *.  P.W  (G..,  -  G,..)  I,,. 


ik 

T 


cos  0(5)  0 

- 2 -  +  - 

dn 


9J.S) 


r^o  Pg(0  sin0/5^)(G„,,  -  G„.i) 

O 


^  2^  cos0g(5)sin0^(5^)(G„,i+G„.i) 

-sin0^(5)cos0g(5^)  GJ 


(159) 


(160) 


56 


Similarly,  the  Fourier  expansions  of  eqs.  (158)  and  (159)  are: 


GJ  4 


ik  d  ri 


4  ds 


«+l 


/  s  5  rl 

/; 


s)  (G.,.-G,.,)  4 


L  J 


4cos0^(5)sin6g(5^)(G„,i+G„.i)]^, 


/C  /  \ 
+  - 


n  d  r  I 


[(^osd  (s)cosd  (sJG„ 


-sin0»sin0j5^)(G„^j  +G„_j)] 


^<T(^)  =  ~2^fo  P  A) 

-  -  —  r  U  p  (5 )  [  A  G  ]  I, 

2  ds  Jo  °  ds  ”  '’" 

<9 

■"  ^  [cOS0/5)cOS0g(5^)G„ 

+  isin0^„(5)sin0/5^)(G„,,  +  G^.^)]  i,„ 


-  +  A  I  r  P  (5 )  (G  j  +  G  j) 

4[  Pg(5)  a«^.Jio  <i>^” 

A7rcos0f5)  ./ 

- ^  r  p  (5 )  (G  ,  -  G  ,)  K, 

4p  (s)  J  o  °  '' 

/I  fcOS0G)  all  r/,  ...... 


4  p<,(.s)  a«5. 


p/sj  smd/sJ(G„^^  -  G„_^) 


/J*.  P.(^.)  [^  COS0/5)  sin 0^(5^)  (G„,i  +  G„_i) 
-sin0^(5)cos0^(5^)  GJ 


(161) 


(162) 
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